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> ' I. INTRODUCTION 

Since the early work of Bochner [1] on the classification of polynomial solutions of second-order linear homogeneous 
differential equations, the problem of finding polynomial solutions to a given differential equation has attracted the 
attention of many researchers Q-j^]. A reason for such interest is that many problems in quantum mechanics, specially 
those arising from Schrodinger's equation, after the separation of the asymptotic-behavior factor of the wave function, 
' there remains a polynomial- type factor in the solution fiol ]- jl8| . Given a differential equation 

y" = X Q {x)y' + s (x)y, (1) 



it may appear to be a simple question to ask whether Eq.([T]) has a polynomial solution or not. In fact, unless the 
given differential equation lies within the classification scheme of Bochner [l[ , the answer to this question seems far 
from simple. The problem seems even much harder if we ask whether Eq.([T]) has polynomial solutions P n (x) not 
necessarily in a single variable x but rather a polynomial- type solution in a function f(x), say 

n 

P n (x)=J2<Xk[f(x)] k - 

k=0 

This problem was recently discussed by Hall et al in the context of finding exact solutions of Schrodinger's equation 
with soft-core Coulomb Potential There are a number of well-known criteria [l[ for determining if there are 

polynomial solutions of the following differential equation 

{a 2 fiX 2 + a 2 ,ix + 02,2)2/' + (aifiX + a M )z/ - T 0fi y = 0, (2) 

which we now outline. 
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Theorem 1. ([14], Theorem 3) The differential equation has a nontrivial polynomial solution of degree n G N 
(the set of nonnegative integers) if, for fixed n, 



T , = n(n - 1) a 2 , + n a 1 ^ 0l 
where the polynomials may be readily obtained from the recurrence relation (n = 0, 1,2. . .): 
((2n + l)a 2 , + ai, )(2(n + l)a 2 ,o + ai,o) 



(3) 



(na 2 ,o + ai,o) 

((2n + l)a 2: o + ai,o)(2n(n + l)a 2:0 a 2 ,i + 2(n + l)ai j0 a 2 4 - 2ai^a 2 ,o + 01,001,1) 

(na 2 ,o + ai,o)(2na 2 , + 01,0) 

(n + l)(2(n + l)a 2 . + ai,o)(4r7, 2 a 2 , 2 a 2 + a^oa^i + 4na 2i oai, a 2 , 2 - n 2 a 2 , a 2 i + a 2 a 2 , 2 - 01,101,00^1 — nai^ a| i) 



2/n+l 



(na 2j0 + ai,o)(2na 2:0 + ai, ) 



(4) 



initiated with 



y = 1, yi = a^oOi + ai,i 



This theorem classifies most of the standard orthogonal polynomials such as Laguerre, Hermite, Legendre, Jacobi, 
Chebyshev (first and second kind), Gegenbauer, and the Hypergeometric type, etc. As a simple illustration of this 
Theorem, we consider the differential equation 



x 3-3 + (2a; + 2)- r ofl y = 0. 

ax ax 



For polynomial solutions, we must have To o = n(n — 1) + 2n = n(n +1), n € N, and the recurrence relation 
generating the polynomial solutions now reads 

2M+2 = 2(2n + 2>)xy n+ x + 4y„, y = l,yi = 2x+ 2. 

These polynomials were studied earlier by Krall and Frink and are known in the literature as Bessel polynomials fl9j | . 

Theorem 2: ([§j], Theorem 3) A polynomial P m {z) — "Y^ 1 j=Q a j z ^ * s a solution of the equation if and only if for 
any z € C (the set of all complex numbers) and n > m + 2 it satisfies the following mean-value formula: 
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Pm{z) 



AnQ 2 {z) + m(m - 1)Q 2 " + 2mQi' ^ 



7 £(Qi(«)A/, + 2Q 2 (z))P m (z + A p 



(5) 



where Q 2 (z) = a 2 ,o z2 + 02, i z + a 2,2, <3i(-z) = ai o^ + 01 1 and £/ie numbers A M = A M (n, 2) are roo£s 0/ i/ie equation 

N 



\ n— 2fi 

y (-1)^^— - = 0, iv 

^ v y 2^m 



(6) 



/or m, fc, and rt oe positive integers such that n > k + m and the square brackets [rj\ denote the integral part of a 
number rj. 

A polynomial P m (x) is a Bessel polynomial if and only if it satisfies the following mean- value formula: 

2((x + l)X^ + x 2 ) 



P m (x) = 



— J 2nx 2 4- m(m + 1) 



P m (x + A M ), n>m + 2 



where A p arc the roots of the equation 



N 



\n—2fi 

£(-iy^-^ = u. v 

fj,=a 
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In the present work, we find the conditions under which the differential equation 

(a 3j0 a; 3 + a 3 .ix 2 + a 3 . 2 x + a 3 , 3 ) y" + (a 2fl x 2 + a 2 .\x + 0,2,2) y' - (n,oX + t^x) y = 0, 



(7) 



has polynomial solutions of the form y n (x) — ^ a kX k . This class of equations contains a number of important 

k=0 

differential equations such as the confluent Hcun equation, biconfiuent Heun equation and the general Heun equation 
[ill-Q, which are usually studied individually in the literature. Some general results concerning solutions of Eq.((7j 
may be found in Ref.[H[. Our method of studying the polynomial solutions of Eq.Q relies on the asymptotic iteration 
method (AIM), which can be summarized by means of the following two theorems. 

Theorem 3: ([13], equations (2.I3)-(2.I4)) Given X = \q{x) and s = so{x) in C°° , the differential equation (QJ) 
has the general solution 



if for some n > 



where, for n > 1, 



y = exp 



C 2 + Ci / exp / (A (r) + 2a(r)) dr dt 



An A„_i 



\i(t)dt 



= a(x), equivalently 5 n (x) = \ n s n -i — A n _is rl = 0, 



A n — A„_j + s„_i + AoA„_i, 
s n = s' n _ 1 + SoA n _l. 



(8) 



(9) 



(10) 



Theorem 4: Theorem 2) (i) If the second order differential equation (QJ) has a polynomial solution of degree n, 

then 

S n (x) = A„S„_i - A„_iS„ = 0, (11) 

where X n and s n are given by I110\) . Conversely (ii) if X n X n -i 5^ and 5 n (x) = 0, then the differential equation (QJj 
has the polynomial solution of degree at most n. 

Theorem 4 may be the simplest criterion available in the literature for testing whether a given differential equation, 
such as |T]), has polynomial solutions or not. 

The conditions given in the present work are considered new as they generalize the earlier work of Krylov and Robnik 
[lll|-[l2l| and they are explicitly written for practical use within the scope of physical applications. They answer 
directly, without further analysis, whether a given differential equation of the form (7) has polynomial solutions or 
not. Our results, especially theorem 5 below, are of important practical use in solving, for example, Schrddinger's 
equation where the construction of the eigenfunctions for bound states, after extracting the asymptotic behaviors, are 
based on the requirement of termination for some hypergeometric series [H|]-[l6j]. Most of the available techniques are 
usually based on constructing classes of differential equations, see for example the construction of Turbiner [f| and 
the work of Krylov and Robnik (l2j ]. that admit polynomial solutions. Very few results are available to test whether 
a given differential equation has actually a polynomial solution. The present results reported in the paper, Theorems 
5 and 6, can be seen as criteria to test the solvability of a differential equation of the form (|7|) in terms of polynomial 
solutions. 

The present article is organized as follows. In the following section we find the conditions under which Eq.([7]) has 
polynomial solutions if there is any. Therein, we discuss some examples found in mathematical physics literature, in 
order to demonstrate the validity of our general approach. In Section (III) we examine the polynomial solutions of 
the confluent, biconfiuent and general Heun equation. In Section (IV) we discuss a special case of a more general class 
of second-order differential equations 



a k+1:i x i \y' - > T kyi x l \y = 0, k = 0, 1, 2, . 





which admit polynomial solutions. 
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II. POLYNOMIAL SOLUTIONS OF EQUATION © 

Theorem 5. The second-order linear differential equation has a polynomial solution of degree n if for any pair 
of the coefficients 03,0; «2,o an d r i,o n °t simultaneously zero, 



n,o = n(n - l)a 3 . Q + na 2 ,o, n = 0, 1, 2, . . . , 
along with the vanishing of (n + 1) x (n + 1)- determinant given by Tabled 
Proof: The proof of this assertion follows by use of theorem 3, equation ([9]), which yields for 



5% = if 7-1,0 = 02.0 and 



n,i -«2,2 

7"1,0 n.l — a 2,l 



= 



(12) 



82 = if n,o = 2ct3.o + 2a 20 and 



7"1,1 —0>2,2 — 2«3,3 

Tl,0 Tl, l~a2 : l — 2<23 : 2 — 2d 2 , 2 

n,o - a 2 , n,i - 2(a 3a + a 2 ,i) 



if no = 6a 3 , + 3a 2 ,o and 



n,i -012,2 -2a 3j3 

n,o n,i - 02,1 -2a 3 , 2 - 2a 2 ,2 -603,3 

n,o - a 2 ,o n,i - 2a 3 ,i - 2a 2 ,i ~6a 3 , 2 - 3a 2 , 2 

n o - 2a 3 , - 2a 2 , n,i - 6a 3 ,i - 3a 2 ,i 



and so on. A procedure which can be easily generalized for S n — to yields n,o = Ti(n — 1)03,0 + ^^2,0 subject to 
vanishing (n + 1) X (n + l)-determinant A„ + i = where A„ + i given by Table |TJ The derivation of this determinant 
can be obtained by substituting f(x) = X)fc=o c k% k to yield the four-term recurrence relation 

[n,o -(k-2)(k- l)a 3 , - (k - l)a2,o] c k -i + [n,i - k((k - l)a 3 ,i + a 2 ,i)] c k -(k+ l)(ka 3 . 2 + a 2 , 2 )c fe+ i 

- (k + 2)(k + l)a 3 , 3 c fc+2 = 0, c_i=0 



from which the determinant follows directly. 

Theorem 6. A necessary condition for the second-order linear differential equation 



'k+2 



fk+l 

E 



„k+l-i 



y 



y = 



\i=0 



to have a polynomial solution of degree n is 

Tfe,o = n(n - l)a k+2 fl + na k+lfll k = 0, 1, 2, . . 



□ 



(13) 



(14) 



Proof: This result follows by putting y(x) = P n (x) where P n {x) = ^2 o,jX J in (|13[) and then multiply out everything 

3=0 

and equating the coefficients of each power of x. The highest power of x yields the necessary condition (TT4")) . O 

Example 1: ( (lTj[ , equation 39) In their discussion of Schrodinger equations that allowing polynomial solutions, 
Krylov and Robnik [llj investigate the solution of the differential equation 



x A y" + a{x z - l)y + {/3x + 7)2/ = 0. 



(15) 



The polynomial solutions of this equation follow directly by use of Theorem 5. Specifically, by using Eq. (IT2l . we find 
the following condition for polynomial solutions 



P = —n 2 — (a — l)n., n = 1, 2, . . . 
where the parameters a, (3 and 7 must satisfy: 



(16) 



TABLE I: The determinant A n +i = for the polynomial solutions of the differential equation (J7J. 

n,i -a 2 ,2 -2a 3 , 3 ... 

71,0 7l,i ^ 02,1 — 2(d3,2 + 122,2) — 6a 3 , 3 ... 

7i,o— <X2,o 7i,i - 2(a 3 ,i + a 2 ,i) — 3(2a 3 , 2 + 02,2) ... 

... 7i,o — (n — 2)(n — 3)a3,o — (n — 2)02,0 n,i — (n — l)(n — 2)03,1 — (n — 1)02,1 — n((n - 1)03,2 + 02,2) 

... ri,o — (n — l)(n — 2)03,0 — (n — 1)02,0 7i,i — n((n — 1)03,1 + 02,1) 



6 



-7 a ... 

-/3 -7 2a ... 

-p - a -7 3a ... 











-/3-(n-l)(n-2)-(n-l)a -7 



= 0. 



In particular, 



n = 1 => /3 = —a and 
n = 2 => P = -2(q + 1) and 



—7 a 
a —7 



7 2 — a 2 = => 7 = ±a, 



—7 a 
-/3 -7 2a 
-0-a -7 



= - 7 (2a 2 + 3a^ + 7 2 ) = 



7 = 0,±V2a(2a + 3). 



Example 2: ([10], Eq. (4)) In their study of the polynomial solutions of the planar Coulomb diamagnetic problem, 
Chhajlany and Malnev [lO] studied the following differential equation 



y" + (p - 2x 2 )y' + (Sx + a)y = 0. 
The equation admits polynomial solution of degree n for 

5 = 2n 

subject to the vanishing of the following determinant 



(17) 
(18) 



-a -p -2 
— S —a —2p —6 
-5 + 2 -a -3p 











-<5 + 2(n-l) -a 



(19) 



Example 3: (15], Eq. (49)) Boztosun et al. studied analytical solutions of the Bohr Hamiltonian for the Davidson 
potential 



(20) 



where /3q is the position of the minimum of the potential (|15||. Eq. 3). Their analysis reduces to the investigation of 
the exact solutions of the differential equation ([l5j, Eq. 49) 



fl L {fi) = (2(3 - ^±H)/; L (/3) + ( 2/ x + 3 - e)f n , L (/3), 

which can re- written, by comparison with Eq. (|T|). as 

(3f"M ~ ' + l ))fn, L {P) - (2/i + 3 - e)f}f n , L {fl) = 0. 
Using Theorem 5, for the polynomial solutions, we must have 

e„ = 2/i + 3 + 4n, n = 0,1,2,.... 
The polynomial solutions using Eq.® are 

f yo(x) = 1 

y x (x) = 2x 2 - 3 - 2^ 

y 2 {x) = 4a; 4 - 4(5 + 2/i).x 2 + (3 + 2 M )(5 + 2/i) 

y 3 (x) = 8x 6 - 12(7 + 2fi)x i + 6(7 + 2//)(5 + 2fi)x 2 - (3 + 2/i)(5 + 2/i)(7 + 2(i) 



(21) 



(22) 
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III. POLYNOMIAL SOLUTIONS OF THE HEUN EQUATION 



A. Confluent Heun's Equation 

The confluent Heun's differential equation [2(|-[24|], written in the simplest uniform shape, is 

P + l 7 + l\ , f fi v 



y" + [a + ■ + 

2 2—1 



y 



y = 



2 2—1 

which can be written as 

z(z - l)y" + (az 2 + (7 + /3 - a + 2)z - a + l) y 1 + ((/i + v)z-fi)y = 0. 



(23) 



(24) 



The polynomial solutions of this differential equation can be found easily using Theorem 5. For polynomial solutions, 
we must have 



/i + v = —na, 

along with the vanishing of the (n + 1) x (n + l)-tridiagonal determinant given by 



(25) 



H a- 1 

na n- (■y + /3-a + 2) 2(a - 1) 

(n-l)a n-2(j + /3 -a + 3) 3(a - 1) 



















a /i — n(n — 1 + 7 + /3 — a + 2) 



(26) 



B. Biconfluent Heun's equation 

For the biconfluent Heun's differential equation [2l[ 



xy" + (~2x 2 -px+(a+ 1)) y' + ((j - a - 2)x - ^(6 + {a + l)0)\ y = 0, 



(27) 



the polynomial solutions can follow again by use of Theorem 5. The condition for polynomial solutions, in this case, 
is explicitly given by 



7 = a + 2(n + l), n = 0,1,2,..., 
along with the vanishing of the (n + 1) x (n + l)-tridiagonal determinant given by Table [TT1 



(28) 



C. Heun's General Equation 



The Heun differential equation [20J- 24| with four regular singular points located at z = 0, 1, a, oo is given by 

'7 S 



V 



y + —1 — 7T7 \y = °> 



z(z — l)(z — a) ' 



(29) 
(30) 



z z — 1 x — a 
where the parameters satisfy the Fuchsian condition 

1 + a + /3 = 7 + 5 + e 

to ensure the regularity of the singular point at infinity. Equation (|29p can be written as 

(z 3 - (1 + a)z 2 + az)y" + ((7 + e + 5)z 2 - (a(<5 + 7) + e + 7)2 + aj) y' + {aj3z - q)y = (31) 

The polynomial solutions of this differential equation follow by use of Theorem 5. For polynomial solutions, we must 
have 



a[3 = —n(n — 1) — n(-y + e + S) 9 J^— a — — n or j3 = — n, 



(32) 



along with the vanishing of the (n + 1) x (n + l)-tridiagonal determinant given explicitly in Table (IIII 



TABLE II: The determinant A n +i for the polynomial solutions of the differential equation (|27(l . 

^{8 + {a + l)/3) -(a + 1) ... 

-( 7 - a- 2) \{5+{a + l)P)+ $ -2(a + 2) ... 

-(7-a-2)+2 i(<5+ (a + !)/?) + 2/3 -3(a + 3) ... 



-(7- a- 2) §(5+(a + l)/3) -n(a + n) 

+2(n-2) -(n-l)(n-2) + (n-l)/3 



_( 7 __ a _ 2 ) (a+l)/3)/3 

+2(n-l) -n(n-l)+n/5 



TABLE III: The determinant A n +i for the polynomial solutions of Heun's equation (|30[1 . 

q -cry ... 

-a/3 q+ (a(8 + j) + £ + 7) -2o(l + 7) ...0 

-q/3 - (7 + e + <5) g + 2(o + 1) + 2(a(5 + 7) + e + 7)) -3(-2(a(<5 + 7) + e + 7) + aj) . . . 

... -a/3- (n- l)[(n-2) g + n(n - 1)(1 + a) 

+( 7 + e + <5)] +n(a(5 + 7) + e + <5) 
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D. Shifted Coulomb Potential 



As an example of a Heun equation that can be treated directly using our results of Theorem 5, we consider the 
Schrodinger equation for a single particle bound by the shifted Coulomb potential 



iP(r) = Ei>(r), Vi(r) 



r + 0' 



f3>0 



(33) 



where A is the d— dimensional Laplacian operator, d > 2. The potential V\(r) has been used as an approximation for 
the potential due to a smeared charge distribution, rather than a point charge, and may be appropriate for describing 
mesonic atoms [2^. This problem has been discussed at length for d = 3 in the two recent articles [13]- [28]. The 
purpose of this section is first to show that the quasi-exact solutions of the eigen-equation (|33[) follows directly from 
Theorem 5, and also to extend the earlier results of Hall et al. (^-[Hl to arbitrary d > 1 dimensions. In atomic 
units, the radial Schrodinger equation (|33|) for the potential V\(r) in d dimensions reads 



2dr 2 



k(k + 1) 



Z 



iP = Ei/j } k = ^(2l + d-3) 



(34) 



2r 2 r + (3 

We may assume the solution of equation Q34[) . which vanishes at the origin and at infinity, is 

where f(r + /3) is to be determined. On substituting Eq. (|3"5|) into Eq. (|3"4")l . we obtain the following second-order 
differential equation for / = f(r + /?) 



(35) 



/" + 2 



(k + 1) 



2a(k + 1) 



2Z 

r- (3 



2E 



/ = o. 



We multiply by r(r — (3) to obtain for E = —a 2 /2 



r[r + 0)f" + [-2ar 2 + 2(k + 1 - afi)r + 2/3(1 + k)] f + [(-2a(k + 1) + 2Z)r - 2a(3(k + 1)] / = 0. 



(36) 



(37) 



Eq. (|37|) is an example of confluent Heun equation [271 ] . however, theorem 5 gives directly the conditions for polynomial- 
type solutions, namely, 



n! =n= 0,1,2, 



fixed 



n' + k + 1 ' 

subject to the conditions on the potential parameters given by the vanishes of the tridiagonal determinant 

72 P% "3 



(38) 







where 



7n-2 (3,i-2 Ctn-1 

7n-l fin-l 0! n 
In Pn 



(3„ = 2a(3(k + n + 1) - n(n + 2k + 1) 
ctn = -nj3(n + 2k + 1) 
7„ = 2a(n - n' — 1) 



(39) 



The meaning of n' is that n' = n but fixed in the sense that for A3 = 0, we should fix n' = 2 for our computations. 
Equation (|3"8"]) gives the exact eigenvalues 



(40) 



subject to conditions on the parameters a, (3 and Z as given by the vanishing of the determinant A„ + i = 0. In Tabic 
TVl we summarize the first few conditions, for n — 1,2,3,4, on the parameter /? in the potential V\(r) = Z/(r + j3) 
so that Eq. (|3"?| has polynomial solution of the form f n (r). It should be noted that theorem 5 tells us explicitly 
whether or not the differential equation has polynomial solutions: finding these polynomials is a problem that in this 
case remains to be solved, for example by using AIM. 
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TABLE IV: Conditions on /3 for the existence of polynomial solutions of Eq. (|37[) , 
a Conditions on j3 
^ °P - 1 = o 

^3 (k + 2)a 2 /3 2 - 3(k + 2)af3 + 2k + 3 = 

(k + 3)(fe + 2)a 3 /3 3 - 6{k + 3)(fc + 2)a 2 /3 2 + (llfc 2 + 50fc + 54)a/3 - 3(fc + 2){2k + 3) = 

jfg (fe + 2)(fe + 3)(fe + 4)a 4 /3 4 - 10(fc + 2)(fc + 3)(fc + 4)a 3 /3 3 + (720 + 823fc + 300fc 2 + 35fc 3 )a 2 /3 2 
-(720 + 925fc + 381fc 2 + 50fc 3 )a,fl + 6(fc + 2)(2fc + 3) (2k + 5) = 



IV. NEW CLASSES OF POLYNOMIAL SOLUTIONS 



The differential equations discussed in the previous sections were characterized by the fact that the parameters Tiq 
and T± t i are the only coefficients in the differential equation that depend on the nonnegative integer n. The parameters 
o-3,i,i = 0,1,2,3 and a,2,j,j = 0,1,2 are constants to be determined, based on the values of T\^,k = 0,1 allowing 
polynomial solutions. In this section we use AIM to discuss new classes of differential equation, where the parameters 
&2,i, i = 0, 1, 2 are also functions of n that allow polynomial solutions. We have 

Theorem 7: For arbitrary values a and b and nonnegative integers of m and n, the solutions of the second-order 
linear differential equation 

„ ax l ~ 2 , m(m + 1)6 
y = + ~T7, ^TT^' * = 2,3, ... (41 



are given by 



™ m+l „ , n ?7i+l 2m + 1 ax 1 1 , . , , 

^ = ^2^! - — —;- , m=l,2,3,...,n = 0,l,2,..., 42 

t — 1 t — 1 I — 1 b(n + m) 



where, for polynomial solutions, n = v(l — 1), for some v = 0,1,2, ... . Here 2 Fi(a, ft; 7; x) refer to the classical Gauss 
hypergeometric series: 

fr[ 7(7 + 1)... (7 + a- i) fci 

In particular, for the differential equation (I = 2): 

11 a 1 m(m + l)b 

y" = j^v 1 + x , {b + ^ ( 44 ) 

the polynomial solutions are given by 

y™=x™+\F 1 (-n,m + l;2(m+l);- °* ), m = 1,2,3, ... ,n = 0,1,2, ... . (45) 

b(m + n) 

Here, each y™ is a polynomial of degree n in the variable x. 

i—i 

Proof: By means of the substitution z = b + a ^ +n , the differential equation Eq.(|4jJ) is reduced to 



d z y 



(n + m) (I- 2) 



(l-l)z (l-l)(z-b) 



dy m(m + 1)6 



dz 2 
Further, let 



dz (l-l) 2 z(z-b) 



-V- (46) 
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we find that the functions f™( z ) now satisfy the differential equation 

f n ( ,_ [{n-m-l)z-{n + m)b] n(m + 1) 

(Z-l)z(z-fe) /l j+ (l-iyz(z-b) J{ h [ <) 

which can easily be solved by use of the asymptotic iteration method, Theorem 3. 

V. CONCLUSION 

This paper is about exact solutions of second-order linear homogeneous differential equations that arise in math- 
ematical physics. If the exact solutions are transcendental functions, then the term 'exact' has a recursive ring to 
it, since such solutions can only be regarded as exact if the solution functions turn out to have been already named 
and studied. An exception occurs when the unknown part of a solution is a factor that is a polynomial, either in the 
independent variable, say x, or in a function of it, f(x). In the study of physical problems where the theory leads to a 
differential equation, such as Schrodinger's equation, it is extremely helpful to have at one's disposal some compactly- 
expressed solutions that typify the kind of exact system behaviour that is implied by the differential equation. It is 
expected that the analysis presented in this paper will provide for a rich and useful variety of such exact solutions. 
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